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Abstract 


The charge density response function and the exchange hole are closely related to each other via the fundamental fluctuation- 
dissipation theorem of physics. A simple approximate model of the static response function is visually compared on several 
examples in order to demonstrate this relationship. This study is completed by illustrating the well-known isomorphism between 
the exchange hole and the square of the dominant localized orbital lying in the space region of the reference point of the exchange 
hole function. The implications of these relationships for the interpretation of common chemical concepts, such as delocalization, 
are discussed. 
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1. Introduction 


To establish links between the wave function and chemical 
concepts characterizing the structure and reactivity of atoms, 
molecules, solids and interfaces has always been a consider¬ 
able challenge for computational and theoretical chemistry. Al¬ 
though this is a subject which has been in the focus of scien¬ 
tists since the very beginning of quantum chemistry, the de¬ 
velopment of Conceptual Quantum Chemistry tools remains 
even in our days a dynamically developing domain, producing 
new ideas and leading to a deeper understanding of old con¬ 
cepts. This research activity betrays a natural need for having 
a sound basis of quantitative or semi-quantitative characteriza¬ 
tion of chemical objects and their interaction in the framework 
of concepts, which are most often deeply rooted in chemical 
thinking. Technically it means that one has to extract specific 
parameters from many-electron wave functions of high com¬ 
plexity, which are able to quantify or illustrate concepts. 

If we are to identify some trends in the evolution of concep¬ 
tual tools over the decades, one of the most remarkable feature 
is that the definition of extracted quantities become more and 
more independent of any specific approximation scheme and 
are rooted more and more in true observables. On the one hand, 
at the beginning, quantities like bond orders and atomic popula¬ 
tions were constructed directly from intermediate quantities of 
popular quantum chemical approximations (molecular orbitals, 
populations analysis in a given basis set, valence bond config¬ 
urations, ...) which could not be interpreted outside the scope 
of a given computational scheme. On the other hand, in the 
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course of the last 20 years, researchers made an effort to define 
generally applicable conceptual measures and indices, which 
are valid independently of any quantum chemical approxima¬ 
tion method. Emblematic examples of such approximation- 
independent concepts are atomic multipoles defined on the ba¬ 
sis of the topological partitioning of the electron density (T), de- 
localization indices between topologically defined atoms (2 - 6| 
or, more recently, the total position spread (TPS) (T), based on 
the localization tensor of Resta CD- Various measures of elec¬ 
tron localization, like the electron localization function (ELF) 
|9 -111, and other analogous functions of the three-dimensional 
space 112 -141 are also defined in a universal manner, without 
referring to specific notions which would be valid only in the 
realm of a particular approximate electronic structure method. 
Such a philosophy has not only the advantage of making pos¬ 
sible a judicious comparison of results obtained from different 
electronic structure methods (e.g. atomic charges obtained from 
Gaussian and plane-wave basis set calculations), but opens the 
way to the comparison of theoretical values directly with the ex¬ 
periment. For instance, atomic charges and Laplacians derived 
from the QTAIM (Quantum Theory of Atoms In Molecules) of 
Bader can be derived not only from computed, but also from 
experimental electron densities obtained from high-resolution 
X-ray diffraction data [ 15 4181. 

Inspired by the works of Parr and his coworkers 03, a 
whole family of functions have been defined in the framework 
of the Conceptual Density Functional Theory, leading to a for¬ 
mally rigorous interpretation of the chemical potential of an 
electronic system as the functional derivative of the electronic 
energy with respect to the number of electrons, and of the linear 
charge density response function as the functional derivative of 
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the electronic energy with respect to the external potential [ 201. 
One can define quantities like hardness, softness, Fukui func¬ 
tion, ... via various higher order and mixed derivatives [2TJ. It 
has to be stressed that these definitions are fully general and 
valid not only in DFT but, in principle, at any level of modern 
electronic structure theory. Among all these quantities, we are 
going to be concerned mainly by the softness, more precisely 
by the softness kernel. 

The relevance of the linear charge density response func¬ 
tion, which is called by some authors linear response kernel 
and which is equal to the negative of the softness kernel, for 
the characterization of localization/delocalization properties of 
an electronic system has been recognized in the past by sev¬ 
eral authors. The remarkable correlation between atom-atom 
softnesses i.e. atomic partition of the linear response functions 
and atom-atom delocalization indices have been pointed out 
through the example of Y-conjugated compounds I22ll23ll . The 
properties of the linear density-density response and its eigen¬ 
value decomposition have been used by Savin and his cowork- 
ers to analyze the strong density-dependence of the Kohn-Sham 
potential in density functional theory lf24l and the radial density- 
density response function has been analyzed for a series of atoms 
along the adiabatic connection path |25l . 

More recently, in a series of papers [ 26-{321, Geerlings and 
his coworkers have analyzed the distribution of the charge den¬ 
sity response in the space for atoms and molecules. One possi¬ 
bility to study this quantity is to calculate the atom-condensed 
linear response parameters, corresponding to the charge-flow 
polarizabilities in the general distributed polarizability theory of 
Stone [33 1 341- Distributed polarizabilities can be implemented 
in the context of various definitions of atoms in molecules (AIM). 
The possibilities cover a wide range of definitions of AIM, go¬ 
ing from the partition of the basis functions to the Hirshfeld- 
like fuzzy <351 or to the QTAIM G) discontinuous partitioning 
of the space. QTAIM charge-flow polarizabilities have been in¬ 
troduced first in |[36l . 

The full non-expanded density-density response kernel, cal¬ 
culated in the uncoupled perturbed Hartree-Fock theory, has 
been plotted in Ref. f30l for a few planar metallic and sim¬ 
ple organic systems in view of characterizing their aromatic¬ 
ity. The connection of the response kernel to the delocalization 
and to induction and resonance effects have been discussed in 
Refs. J2glHI. 

In the present contribution, we discuss not only a few more 
examples to illustrate the behavior of the linear density-density 
response kernel for molecular systems as plots of the charge 
deformation due to a perturbation in a fixed point in the space, 
but we attempt to point out and to demonstrate the possible re¬ 
lationships between the response kernel and the exchange (or 
Fermi) hole function. Furthermore, the well-known connection 
between the exchange hole and localized orbitals is also illus¬ 
trated. 


2. Basic relationships 

2.7. Exact response and exchange-correlation hole functions 

Using the definition of the two-particle (pair) density oper¬ 
ator # 2 ( 11 , 12 ) i n terms of the one-particle density operator h(x_ x ) 
and of d(r 1? r 2 ), the Dirac delta function defined later: 

#2(u,r 2 ) = wfe)wfe) - dfe,^)#^), (1) 

and the definition of the exchange-correlation hole in terms of 
the expectation values n(r_\) and n 2 (r 1? r 2 ): 

h xc (l i,r 2 ) = ——- n( r 2 ), (2) 

*(l 1 ) 

one can deduce the following relationship: 

n (Zi)h xc (li,Z2) = * 2 ( 11 , 12 ) ~ *(li)*(h) 

= (* 2 ( 11 , 12 )) ~ (*(h)) 

= <#(ri)#(r 2 )> - S(r l9 r 2 ) (fife)) - (fife)) (#(r 2 )) 
= <d#(ri)d#(r 2 )> - d(r 1? r 2 ) (fife)) , (3) 

where (O) = ( v Fol O I'Fo) denotes a ground state expectation 
value, and we use the decomposition of the density operator as 
a fluctuation 5h(x x ) around its mean: #(rQ = n(x x ) + Sh( rQ. The 
Dirac delta function S( r l9 r 2 ) = S(r 2 - rQ is usually defined in 
the most general terms via its sifting property f37l 

f dr 1 f(i 2 )S(T 2 -r 1 )=f(r 1 ). (4) 

There are numerous options to represent S(r 2 - rQ by the 
limiting value of a series of functions, like in Eq. ( [13] ). 

Eq. ([3]) leads to: 

(8h(r_ x )8h(x_ 2 )) = 5(x_ x ,x_ 2 )n(x_ x ) + n(x_ x )h xc (x_ x ,x_ 2 ). (5) 

The right-hand side of the above equation is usually referred 
to as the static form factor (cf. Ref. (38)). As we can see, 
Eq. ([5]) establishes a relationship between the density-weighted 
hole function and the fluctuations of the charge density, fluctu¬ 
ations that are themselves related to the response function by 
the virtue of the (zero-temperature) fluctuation-dissipation the¬ 
orem. In the following we derive an approximate relationship 
between Eq. 0 and the static response function. 

We proceed by the technique of effective denominators, usu¬ 
ally related to the name of Unsold |[39l and which has been later 
generalized by others, e.g. MM- The aim of this technique 
is to replace the state-specific energy denominator in the sum- 
over-states expression of the response function by a constant 
and to transform the sum over the excited states in the numera¬ 
tor to a ground state expectation value by invoking the resolu¬ 
tion of identity. 

Using the shorthand notation for the transition density as¬ 
sociated with the a excited state n a (r_\) = ('Fol *(l i) Y^a) and 
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Figure 1: Sectional drawings, for the molecules of ethylene (top left), butadiene (top right), water (bottom left) and naphthalene (bottom right), along three 
orthogonal planes of the isocontours of (1) the response function corresponding to the reference point r 0 indicated by the black dot labelled "Q" and of (2) the 
mathematical object defined by the right-hand side of equation Eq. {8|, shown for the same reference point r 0 . The isocontour levels are in the range ±0.002 for the 
response function and static form factor of the ethylene and butadiene molecules, in the range ±0.005 for the response function and static form factor of the water 
molcule, in the range ±0.0005 for the response function of the naphthalene and ±0.00025 for the static form factor of the naphthalene. The isocontours are separated 
by 0.0005 in the cases of the ethylene, butadiene and water molecule, by 0.0001 for the response function of the naphthalene and 0.00005 for its static form factor. 
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with the excitation energy aj a , the static response function can The summation over the virtual orbitals is transformed into 

be written as: a sum over all the orbitals minus a sum over occupied orbitals: 


Vli,r 2 ;0) = 2^ 


w«(li)w«(r 2 ) 


a^O 

2 


(D n 




^n a (ri)n a (r 2 ) 


a± 0 


(Sh(r t )Sh(r 2 )) 


"(£i,r 2 ) 

2 (Sh(T t )Sh(l 2 )), 


CD 


( 6 ) 

(7) 


where, disregarded the hypothesis of working with real wave 
functions, only the last step involves an approximation. While 
the position-dependent effective denominator function (d(x_\ , r 2 ) 
in line ([6]) is able to maintain the equality, its reduction to a 
position-independent constant co necessarily leads to an approx¬ 
imate expression of the response function. We have used simi¬ 
lar approximations in earlier works [23. 381. Comparing equa¬ 
tions Eq. 0 and Eq. 0 allows us to write: 


occ all 

"(£i,£ 2 )*o(£i,£ 2 ) = 2 ZZ 0p(£ 2 M'(£ 2 M' (£1 ) 0 /j(£i ) 

i p 
occ occ 

-*ZZ 

i j 

where fj is an occupied orbital and designates an arbitrary 
(virtual or occupied) orbital. In this last expression, we rec¬ 
ognize the usual definition of the exchange hole h x (r l9 T 2 ) for a 
single determinant wave function: 


ra(:i)Mli,r 2 ) = -2E^*(£iWEi)^}(£2M(r 2 )> (12) 

ij 

and the resolution-of-identity expression of the Dirac delta func¬ 
tion: 


CD 

^(£i>£ 2 ;0) = £(r|,r 2 )n(ri) + n(r 1 )h xc (r 1 ,T 2 ), (8) 

which is the desired (approximate) relationship between the 
static response function and the right-hand side of Eq. ([5]), the 
static form factor. Note that the two contributions to the static 
form factor, corresponding to the density multiplied, on the one 
hand, by the Dirac delta function and on the other hand, by the 
xc-hole function, ensure the correct charge conservation sum- 
rule. Indeed, the response function on the left-hand side of 
Eq. ([8]) integrates to zero, while, on the right hand side, the 
Dirac delta function integrates to plus one and the exchange- 
correlation hole function to minus one. 

In the following we are going to simplify the problem and 
instead of working with the exact response function and the 
exchange-correlation hole function, we are going to consider 
the noninteracting (or bare) response function and the exchange 
hole function, which arises in an independent particle model. 

2.2. Noninteracting response and exchange hole functions 
In this subsection we are going to apply the effective de¬ 
nominator technique to the zero-frequency bare or noninteract¬ 
ing response function: 

, , ^ C(£ 2 M(£ 2 )€(£iM.(£i) 

xoih’ii) = 2 2_t 2 j -;rr 7 -> < 9 > 

i a a i 

where 0;, e t (0 a , e a ) are occupied (virtual) orbitals and orbital 
energies. Using a position-dependent effective denominator, 
6L>(ri,r 2 ), this yields: 

occ virt 

^(£i>£ 2 )/to(ri,I 2 ) = 2 zz 0a(£ 2 M(£ 2 )0?(£ l)0a(£l)- (10) 

i a 


<1(£i>£ 2 ) = (13) 

P 

Note that the above representation of the Dirac delta func¬ 
tion holds strictly only in the case of a complete orbital basis. 
We obtain: 

"(£i,£ 2 )*o(£i,£ 2 ) = 2S(r t , r 2 ) ^ <P,(r 2 )<P* (£,) + n(r, )h x (r t , r 2 ), 

(14) 

where we can use the Dirac delta function to see the expression 
of the charge density n(x x ) emerge in the first term: 

u(h,l 2 )Xo(h>l 2 ) = £(li,r 2 Mri) + n(x_ x )h x (x_ l9 x_ 2 ). (15) 

As seen before, the simplest approximation to this poten¬ 
tially exact relation is to replace oj( x x , r 2 ) by a position-inde¬ 
pendent constant, Td\ 

^To(li,r 2 ) * d(r 1? r 2 K x x ) + nix^h^x^). (16) 

Compared to Eq. 0. which relates the exact interacting 
static response function to the exchange-correlation hole, equa¬ 
tion ( fl6| ) involves the noninteracting static response function 
and the exchange hole. 

Our average energy denominator approximation is expected 
to fail in reproducing the nodal structure of the response func¬ 
tion in all its details, but the major features are expected to be 
preserved. In the present article all the numerical examples will 
be limited to the noninteracting (bare) response function and to 
the corresponding exchange hole. 
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2.3. Exchange-hole and localized orbitals 

The idea that localized orbitals are closely related to the ex¬ 
change hole has probably first appeared in the work of Luken 
[43], 44) at the beginning of the eighties and it has been later re¬ 
iterated by Tschinke and Ziegler E) as well as by others (46). 
The central idea is that the definition of the exchange hole for 
a single determinant wave function of Eq. ( [12] ) can be consid¬ 
erably simplified if it were possible to use strictly localized or¬ 
bitals. Indeed, strictly localized orbitals are such that the "dif¬ 
ferential overlap" 0* (h)<f>j(Zi) is (almost) zero for i ± j and it is 
|0*(r 1 )| 2 for i = j. Therefore we can write the hole function as: 

Mr 1 ,r 2 ) = -2j]^|i|0,(r 2 )| 2 . (17) 

Let us consider a partition of the physical space to approx¬ 
imately disjoint regions Q; ED, each belonging to a strictly 
localized orbital 0* . In every region only one localized orbital 
contributes to the electron density, i.e. ^(r^l ~ n(r i). Using 
the window function 0/(rj) (which is equal to unity if h is in 
£li and zero otherwise) it is possible to write that: 

hx(l\,h) = -2 Z k‘fe)| 2 ’ ( 18 > 

i 

i.e. that the exchange hole associated with a reference point r x 
is approximated as the negative of the square of the localized 
orbital k which is situated in the space domain of this point and 
satisfies ®k(li) = 1, that is to say: 

Wli>I 2 ) = —2|0fe(r 2 )| 2 forr 1; r 2 e D. k . (19) 

In Section [3j it is shown in several illustrations how far this 
reasoning holds in different systems characterized by different 
degree of localization of their electrons. 

3. Picturial comparison of response functions, 
exchange holes and localized orbitals 

3.1. Computational details 

A software has been written to calculate the two-variable 
functions Xo(lk> foX £(&> lo) and h x (r k , r 0 ) appearing in equa¬ 
tions and ( [12] ) on regular grids {r k }, for a fixed ref¬ 

erence point r 0 , permitting us to generate visual representations 
of both sides of Eq. ^ and Eq. ( fl9| ). The program takes as 
input the orbitals and the total electronic density of a system 
in the CUBE file format using the MOLPRO program (47) . The 
calculations have been done at the Kohn-Sham level using the 
LDA functional and the aug-cc-aVTZ basis set. The sum rules 
f dr k xoilk’ lo) = 0 , / dx_ k h x (r k ,x M ) = -1 and f dv k <5(r*,r 0 ) = 1 
have been checked to verify the suitability of the grid by sim¬ 
ple summation over the grid points. The Dirac delta function 
has been represented by a spherical Gaussian model: 6{r k , r 0 ) « 
^. e (-i w2 \lk~lo\ 2 )' The half-width parameter w = 2.5 has been 
found on our grids to produce a normalization integral closest 
to 1 for all studied systems. Localized orbitals were generated 
using the Foster-Boys localization criterium [48, 49]. 


3.2. Response function and exchange hole 

In order to visualize the relationship between the response 
function and the hole function, we present contour plots of the 
noninteracting response function, 

To (&, r 0 ; 0) and 6(r k , r 0 )n(r k ) + n(r k )h x (r k , r 0 ), 

in cross sections of the molecule lying in the xy, yz and zx 
planes. 

Figure [I] shows the correspondance between the static re¬ 
sponse function x(lk ,r 0 ; 0) and the static form factor, i.e. the 
function appearing on the right-hand side of Eq. 0- We have 
used as example the ethylene, butadiene, water and naphtha¬ 
lene molecules. The reference points r 0 are labelled with a "Q" 
on the Figures. The quality of the ressemblance between the 
static response function and the static form factor does not de¬ 
pend on the choice of the reference points, but these are placed 
on symmetry planes or axis for each molecule and are chosen 
such that they coincide with a grid point. Due to the relatively 
rough grids, an arbitrarily placed reference point would have 
destroyed the symmetry of the resulting figure. For the ethy¬ 
lene and butadiene molecule the reference point is near the C-C 
single bond, outside of the plane of the molecule. In the case of 
the water molecule, the reference point is on the bisector of the 
H-O-H angle. The reference point for the naphthalene molecule 
is near the central C-C bond, above the molecular plane. 

Overall, the correspondance between the two functions is 
of rather acceptable visual quality, in spite of some discrepan¬ 
cies, which can be attributed to the relative crude model of the 
response function, supposing a simple, position independent, 
proportionality with the form factor. The incompleteness of 
the model can be conjectured from the fact that linear response 
functions possess a richer topological structure, as compared 
to the form factor based model. The mathematical structure of 
the form factor makes clear that the dominant feature in both 
functions should be a positive peak centered on the reference 
point r 0 and an essentially negative region corresponding to the 
exchange hole. 

The comparison of the butadiene and the naphthalene clearly 
shows the difference in the extent of the delocalization in these 
two systems. While in the former case the response function 
describes an induced change of the charge which extends out 
only slightly from a central perturbation to the terminal car¬ 
bon atoms, in the case of the naphthalene the charge density is 
quite significantly perturbed on the entire molecular skeleton. 
We can observe also an alternating sign of the charge density 
variations. The topology of the corresponding independent par¬ 
ticle form factor function shows a reasonable resemblance with 
the behavior of the linear response function, although some fea¬ 
tures seem to be different on the two contour plots. For instance, 
the multiple maxima observed in the response function plot of 
the naphthalene molecule seems to be absent on the form factor 
plot. In the case of the butadiene the strong positive peaks at the 
on the opposite side of the molecular plane with respect to the 
reference point under the Cl and C2 carbon atoms are almost 
invisible on the form factor plot. 
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Figure 2: Sectional drawings, for the molecules of ethylene (top left), butadiene (top right), water (bottom left) and naphthalene (bottom right), along three 
orthogonal planes of the isocontours of (1) the exchange hole corresponding to the reference point r 0 indicated by the black dot labelled "Q" and of (2) the square 
of the localized orbital which centroid is situated in this space domain. The isocontours are in the range ±0.03 and separated by 0.005 for the ethylene molecule, in 
the range +0.05 and separated by 0.005 for the butadiene, and in the range +0.005 and separated by 0.0005 for both the water and naphthalene molecule. 
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Figure 3: Effect of the choice of the reference point, ro, denoted by Q4 and Q7, respectively for the water (left) and ethylene (right) molecules. The reference point 
is moving outwards along a symmetry axis of the molecules, remaining in the same bond domain. Column (1) represents the response function in the xy plane, 
column (2) illustrate the mathematical object defined by the rhs of Eq. {8j, both shown for the reference point ro- 


The contour plots in Figure [T] depend on the choice of the 
reference point. We have examined the sensitivity of the re¬ 
sponse function and of the form factor with respect to the posi¬ 
tion of the reference point. Since the shape of these functions 
is expected to change radically if the reference point is moved 
from one electron pair domain to another, we focused our atten¬ 
tion to relatively small displacements within the same domain. 
As illustrated by Figure [3j the shape of the response function 
does not change appreciably; we have the positive peak, associ¬ 
ated in the form factor model with the Dirac delta function and 
a negative region related to the exchange hole in the form factor 
model. However, according to the selected examples, the water 
and ethylene molecules, some details may change in the case of 
a shift of the reference point along the symmetry axis and ad¬ 
ditional critical points may appear in the negative (red) regions 
of both the response function and the form factor model. 

3.3. Exchange hole and localized orbitals 

Figure [2] illustrates the relationship seen Eq. ( [T9] ) between 
the localized orbital and the exchange hole, by plotting 

hxih’lo) and ~\Mlk)\ 2 - 

The reference points r 0 for the exchange hole have been 
chosen at the centroids of the localized orbital obtained 
by the Foster-Boys localization criterium. The following or¬ 
bitals were selected: the C-C bonds for the case of the ethylene 
molecule, the C-C simple bond of the butadiene molecule, one 


of the oxygen lone pairs for the water molecule and the cen¬ 
tral C-C bond of the naphthalene molecule. The grids are con¬ 
structed such that the reference points coincide with one of the 
grid points. 

In the case of the small water and ethylene molecules and 
even in butadiene the contour plots of the exchange hole and of 
the square of the localized orbital are practically indistinguish¬ 
able. Although the agreement is very good in naphthalene as 
well, the squared localized orbital seems to be slightly more 
extended on the naphthalene rings than the exchange hole. 

It has been observed for a long time that bond polarizabil¬ 
ities are remarkably transferable from one system to another. 
This point was discussed from the theoretical chemistry view¬ 
point by Claverie [50] by identifying bond polarizabilities and 
localized orbital polarizabilities. The isomorphism of the local¬ 
ized orbitals and of the exchange hole, which in turn is related to 
the local charge density response puts this observation in a new 
context a provides a sort of explanation in terms of the model 
independent exchange hole through its close relationship with 
the linear response function itself. 

3.4. Discussion 

From a conceptual point of view it is quite obvious that the 
notion of delocalization is primarily connected to the question 
"How electrons respond to an external disturbance?". The ba¬ 
sic physical observable in this respect is the charge density re¬ 
sponse function, unfortunately its direct use in routine studies 
is made difficult by some obstacles. 
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The main difficulty is that the response function depends 
on two space variables so its full pictorial representation in 
the three dimensional space or in any lower dimensional cross- 
section (plane or line) should be done as a function of the ref¬ 
erence point, which should, in principle, scan the whole space. 
Geerlings and his co-workers have recently published a series 
of articles, which present and discuss response functions (lin¬ 
ear response kernels) for a series of systems 132]. These studies 
have been done most often on the noninteracting response func¬ 
tion, which requires significantly less computational resources, 
than the analogous plots of the interacting response, e.g. at the 
coupled perturbed Kohn-Sham or Hartree-Fock or correlated 
(coupled cluster) level. 

A way to get rid of the second variable would be to construct 
domain-averaged quantities, in analogy to the domain-averaged 
Fermi holes (DAFH) EMI : 

Xn A (li)= f * 2 ^( 11 ,r 2 ). (20) 

Jci A 

To the best of our knowledge this quantity constructed from 
the response function has not yet been proposed in the liter¬ 
ature. For instance, it is expected that the shape of the atomic 
domain-averaged response function provides information about 
space regions where the electrons of this particular atom are 
expected to be delocalized. A detailed study of such domain- 
averaged response functions will be be the subject of a forth¬ 
coming publication. 

A further simplification in the characterization of the linear 
response function can be achieved by a double domain average, 
leading to a quantity identical to the charge-flow polarizability 
in the distributed polarizability formalism: 

Xn A ,n B = dv_ 2 I (21) 

To a ,o b measures the propensity of the charge density in £l A 
to flow towards the domain Q B under the effect of an electric po¬ 
tential difference between the two domains. Application of the 
above definition to QTAIM domains defined by Bader (T) leads 
to atom-atom charge flow polarizabilities [36, EH- Although 
relatively small systems were studied in these works, clear sig¬ 
natures of strong delocalization have been observed for the di¬ 
cyan (NCCN) and benzene molecules. In the former case, the 
end-to-end charge flow is remarkably high (-0.499) which is to 
be compared to the charge flow polarizability between bonded 
nitrogen and carbon atoms (-1.256) while the charge flow be¬ 
tween non bonded carbon and nitrogen atoms is very small (- 
0.042). In the benzene, the strong para charge-flow is remark¬ 
able (-0.316) and the opposite sign meta charge flow (+0.103). 
Such a behavior is in good agreement with well-known reactiv¬ 
ity rules in substituted benzenes. 

Our attempt to model the main features of the response 
function in terms of the density-weighted exchange hole was 
motivated in a great extent by the observation that the charge 
flow derived from the linear response function by a double atomic 
domain average correlates well with delocalization indices, which 
can be calculated from the exchange hole. Although the model 


presented in this work does correctly reflect the main trends, it 
can and should be improved in the future. One possible way 
consists in finding a simple model for the position-dependent 
effective excitation energies. Work in this direction is in progress. 

A recent exhaustive review by Geerlings and his coworkers 
l32l came to our knowledge after completing our work. Their 
paper discusses in depth various chemical applications of the 
linear response function, like mesomeric effects, delocalization, 
aromaticity, reactivity, etc of the response functions and quan¬ 
tities derived therefrom. 

One of the roles of the density-weighted exchange hole is to 
establish a link between the most general delocalization mea¬ 
sures, the response function and the localized orbitals, which 
are usually considered as outdated. The quality of the conven¬ 
tional localized orbitals depends on the localization criteria. It 
is interesting to notice that one can derive one-electron func¬ 
tions, which play an analogous role as the localized orbitals, di¬ 
rectly from the exchange hole. In this respect, one should men¬ 
tion the Fermi hole based localized orbitals of Luken and Cul¬ 
berson El El, the natural orbitals derived from the domain- 
averaged Fermi hole by Ponec (53) . the domain natural orbitals 
(DNO) 1521, which were further generalized by diagonalization 
of the ft-th order cumulant density matrices, leading to Natu¬ 
ral Adaptive Orbitals (NAdOs) ED. In this respect one should 
mention Cioslowski’s isopycnic transformation method to ob¬ 
tain localized natural orbitals from correlated (many-determi- 
nantal) wave functions, based on the invariance of the Fermi 
hole [[58]. 

4. Conclusions 

The basic physical model that the (de)localization of elec¬ 
trons is intimately related to fundamental physical observables 
like the charge density response function of the system on the 
one hand and to the exchange-correlation hole function on the 
other, has been illustrated on a few selected examples. Note that 
in independent particle theories we can consider the exchange 
hole and the noninteracting linear response function, which are 
much simpler to calculate than the full linear response and the 
full many-body exchange-correlation hole. 

Pictures of two-dimensional cross-sections of the molecu¬ 
lar space provide a detailed insight to the relationships which 
have been established earlier between charge-flow polarizabil¬ 
ities and atom-atom delocalization indices. Further and more 
rigorous comparisons should be done in the future, since the 
visual shape of the contour plots may be strongly dependent 
on the choice of the contours, selected here on a linear scale 
and uniformly for all the presented pictures. The still objec¬ 
tively existing discrepancies are expected to be removed by an 
improved, position-dependent effective excitation energy mod¬ 
ulating the form factor function, as indicated previously. 

The relationship between localized orbitals and the exchange 
hole, first observed by Luken more than 30 years ago, permits 
an a posteriori justification of the use of localized orbitals in 
qualitative interpretations of bonding and localization. In the 
mean time this analysis underlines some inherent limitations of 
the localized orbital picture. 
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There are several physical implications of the relationships 
which we attempted to make more plausible via graphical illus¬ 
trations. First, we can see that the charge density of an elec¬ 
tronic system at point r 2 responding to an external perturbation 
applied in r_\ can be quite reasonably predicted from the ex¬ 
change hole. Roughly speaking, the response will be nonzero 
essentially in those region where the exchange hole is nonzero 
too. Furthermore, it may be surprising that the charge density 
response in genuinely localized systems remains local in real 
space: it extends only to a few-atom region around the per¬ 
turbation. An implication of this observation is that the charge- 
flow between such regions, which are in close resemblance with 
the exchange hole and by consequence to localized orbital do¬ 
mains, can be relatively small. Selecting such domains for a 
multi-center, distributed description of intermolecular forces, in 
particular in the case of induction and dispersion interactions, 
leads naturally to models without significant charge-flow con¬ 
tributions. In contrast, in metal-like strongly delocalized sys¬ 
tems the possible necessity of including charge-flow contribu¬ 
tions lf59l becomes obvious by the failure of finding localized 
orbitals or in more general terms, by the inherently delocal¬ 
ized nature of the density-weighted exchange hole. Since the 
Resta localization tensor 00 or using a recently suggested al¬ 
ternative name 0, the Total Position Spread Tensor (TPST), 
is the second moment of the density-weighted exchange hole, 
one can establish also a direct link |[60l to the well-known near¬ 
sightedness concept of Walter Kohn |6F, 6ll, which provides a 
general framework to discuss localization and delocalization in 
solids but also in finite molecular systems. 
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